1 Introduction

This paper contains estimates for the effective reproduction number \(R_{t,m}\) over time \(t\) in various provinces (and health regions) \(m\) of Canada. This is done using the methodology as described in [1]. These have been implemented in R using EpiEstim package [2] which is what is used here. The methodology and assumptions are described in more detail here.

This paper and it’s results should be updated roughly daily and is available online.

As this paper is updated over time this section will summarise significant changes. The code producing this paper is tracked using Git. The Git commit hash for this project at the time of generating this paper was 3a5554873d47a90556cbb9b5bad9159dce6322d0.

2 Data

Data is downloaded from the Git repository associated with [3]. This contains the daily cases and deaths reported for Canada by province and health regions.

Repatriated cases are removed and provinces as found in [3] are grouped per the table below:

Province Groupings
Province Province Grouping
Alberta Alberta
British Columbia British Columbia
Manitoba Manitoba
New Brunswick Maritime Provinces
Newfoundland and Labrador Maritime Provinces
Nova Scotia Maritime Provinces
Nunavut Territories
Northwest Territories Territories
Ontario Ontario
Prince Edward Island Maritime Provinces
Quebec Quebec
Saskatchewan Saskatchewan
Yukon Territories

The only adjustments to the data relates to large numbers of deaths reported in Ontario on 2 and 3 October 2020 “due to a data review”. 111 of the deaths reported on these two dates are deaths that occurred during the prior spring and summer (see [4] and [5]). Based on this these deaths were removed from 2 and 3 October and added back in prior days in proportion to deaths reported up un till 30 September 2020. The net effect is no change in reported deaths, but a peak in October is avoided which would have biased estimates.

3 Methodology

The methodology is described in detail here.

4 Results

4.1 Cases and Deaths

Below we plot cumulative case count on a log scale by province:

Below we plot the cumulative deaths by province on a log scale:

4.2 Current \(R_{t,m}\) estimates by Province

Below current (last weekly) \(R_{t,m}\) estimates are tabulated.

Estimated Effective Reproduction Number by Province
province Estimated Type Count (Last Week) Week Ending R - Lower CI R - Mean R - Uppper CI
Alberta cases 2,408 2021-03-05 0.9 1.0 1.0
Alberta deaths 36 2021-03-05 0.5 0.7 0.9
British Columbia cases 3,913 2021-03-05 1.1 1.1 1.2
British Columbia deaths 25 2021-03-05 0.6 0.9 1.3
Manitoba cases 385 2021-03-05 0.7 0.8 0.9
Manitoba deaths 15 2021-03-05 0.8 1.4 2.2
Maritime Provinces cases 81 2021-03-05 0.5 0.6 0.8
Maritime Provinces deaths 3 2021-03-05 0.3 1.1 2.4
Territories cases 17 2021-03-05 0.7 1.1 1.8
Ontario cases 7,406 2021-03-05 0.9 0.9 1.0
Ontario deaths 98 2021-03-05 0.7 0.8 1.0
Quebec cases 5,030 2021-03-05 0.9 0.9 1.0
Quebec deaths 83 2021-03-05 0.8 1.0 1.2
Saskatchewan cases 1,088 2021-03-05 1.0 1.0 1.1
Saskatchewan deaths 13 2021-03-05 0.5 1.0 1.6
Canada cases 20,328 2021-03-05 1.0 1.0 1.0
Canada deaths 273 2021-03-05 0.8 0.9 1.0
Estimated Effective Reproduction Number by Province

Estimated Effective Reproduction Number by Province

4.3 Estimated Effective Reproduction Number for Canada over Time

Below we plot results for Canada as a whole.

Estimated Effective Reproduction Number for Canada over Time

Estimated Effective Reproduction Number for Canada over Time

4.4 Estimated Effective Reproduction Number for Provinces over Time

Below we plot results for each province. We filter out weeks where the upper end of confidence interval for \(R_{t,m}\) exceeds 4.

4.4.1 Alberta

4.4.2 British Columbia

4.4.3 Manitoba

4.4.4 Maritime Provinces

4.4.5 Territories

4.4.6 Ontario

4.4.7 Quebec

4.4.8 Saskatchewan

4.4.9 Canada

4.5 Detailed Results

Detailed output for all provinces are saved to a comma-separated value file. The file can be found here.

5 Discussion

Limitation of this method to estimate \(R_{t,m}\) are noted in [1]

  • It’s sensitive to changes in transmissibility, changes in contact patterns, depletion of the susceptible population and control measures.
  • It relies on an assumed generation interval assumptions.
  • The size of the time window can affect the volatility of results.
  • Results are time lagged with regards to true infection, more so in the case of the use of deaths.
  • It’s sensitive to changes in case (or death) detection.
  • The generation interval may change over time.

Further to the above the estimates are made under assumption that the cases and deaths are reported consistently over time. For cases this means that testing needs to be at similar levels and reported with similar lag. Should these change rapidly over an interval of a few weeks the above estimates of the effective reproduction numbers would be biased. For example a rapid expansion of testing over the last 3 weeks would results in overestimating recent effective reproduction numbers. Similarly any changes in reporting (over time and underreporting) of deaths would also bias estimates of the reproduction number estimated using deaths.

Estimates for the reproduction number are plotted in time period in which the relevant measure is recorded. Though in reality the infections giving rise to those estimates would have occurred roughly between a week to 4 weeks earlier depending on whether it was cases or deaths. These figures have not been shifted back.

Despite these limitation we believe the ease of calculation of this method and the ability to use multiple sources makes it useful as a monitoring tool.

6 Author

This report was prepared by Louis Rossouw. Please get in contact with Louis Rossouw if you have comments or wish to receive this regularly.

Louis Rossouw
Head of Research & Analytics
Gen Re | Life/Health Canada, South Africa, Australia, NZ, UK & Ireland
Email: LRossouw@GenRe.com Mobile: +27 71 355 2550

The views in this document represents that of the author and may not represent those of Gen Re. Also note that given the significant uncertainty involved with the parameters, data and methodology care should be taken with these numbers and any use of these numbers.

References

[1] A. Cori, N. M. Ferguson, C. Fraser, and S. Cauchemez, “A new framework and software to estimate time-varying reproduction numbers during epidemics,” American Journal of Epidemiology, vol. 178, no. 9, pp. 1505–1512, Sep. 2013, doi: 10.1093/aje/kwt133. [Online]. Available: https://doi.org/10.1093/aje/kwt133

[2] A. Cori, EpiEstim: A package to estimate time varying reproduction numbers from epidemic curves. 2013 [Online]. Available: https://CRAN.R-project.org/package=EpiEstim

[3] I. Berry, J.-P. R. Soucy, A. Tuite, and D. Fisman, “Open access epidemiologic data and an interactive dashboard to monitor the COVID-19 outbreak in Canada,” Canadian Medical Association Journal, vol. 192, no. 15, pp. E420–E420, Apr. 2020, doi: 10.1503/cmaj.75262. [Online]. Available: https://www.cmaj.ca/content/192/15/E420

[4] G. Rodrigues, “Ontario reports new record of 732 coronavirus cases, adds 76 more deaths due to data cleanup.” [Online]. Available: https://globalnews.ca/news/7373691/ontario-coronavirus-cases-october-2-covid19/. [Accessed: 01-Nov-2020]

[5] R. Rocca, “Ontario reports 653 coronavirus cases after record number of tests completed.” [Online]. Available: https://globalnews.ca/news/7376209/ontario-coronavirus-cases-oct-3-covid19/. [Accessed: 01-Nov-2020]